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ABSTRACT 

Thermal mid-infrared emission of quasars requires an obscuring structure that 
can be modeled as a magneto-hydrodynamic wind in which radiation pressure on 
dust shapes the outflow. We have taken the dusty wind models presented by Keating 
and collaborators that generated quasar mid-infrared spectral energy distributions 
(SEDs), and explored their properties (such as geometry, opening angle, and ionic 
column densities) as a function of Eddington ratio and X-ray weakness. In addition, 
we present new models with a range of magnetic field strengths and column densities 
of the dust-free shielding gas interior to the dusty wind. We find this family of models - 
with input parameters tuned to accurately match the observed mid-IR power in quasar 
SEDs - provides reasonable values of the Type 1 fraction of quasars and the column 
densities of warm absorber gas, though it does not explain a purely luminosity- 
dependent covering fraction for either. Eurthermore, we provide predictions of the 
cumulative distribution oi E(B — V) values of quasars from extinction by the wind 
and the shape of the wind as imaged in the mid-infrared. Within the framework of this 
model, we predict that the strength of the near-infrared bump from hot dust emission 
will be correlated primarily with L /LEdd rather than luminosity alone, with scatter 
induced by the distribution of magnetic field strengths. The empirical successes and 
shortcomings of these models warrant further investigations into the composition and 
behaviour of dust and the nature of magnetic fields in the vicinity of actively accreting 
supermassive black holes. 

Key words: quasars: general, galaxies: active, galaxies: Seyfert, (magnetohydrody¬ 
namics) MHD, infrared: galaxies 


1 INTRODUCTION 

Researchers currently studying active galactic nuclei agree 
on a few key premises. The central engine is powered by a 
supermassive black hole, and the optical through ultraviolet 
continuum emission is generated by optically thick material 
in an accretion disk within the central parsec down to a few 
gravitational radii. Important specifics such as the thickness 
of the disk as a function of radius are still being resolved, 
but the overall picture is consistent; it is hard to avoid 
disk formation given the persistence of angular momentum. 
Within the central parsec, ionized gas moving at high ve¬ 
locities covers approximately ~ 10% of the sky and gives 
rise to broad (1000s km s“^) resonance and semi-forbidden 
emission lines that characterize Type 1 quasars and Seyfert 
galaxies. In many luminous quasars, the broad UV emission 
lines show signatures of winds; the most obvious of these 
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are the P Cygni-type profiles seen in the high ionization 
reson ance lines of Broad A bsorption Line (BAL) quasars 
fe.g.. IWevmann et al.l[l99lll . More subtle are the blue-shifts 
and asymmet ries seen in the sam e lines in other luminous 
quasars (e.g., [Richards et al.l[2011^ . The wind manifested in 
these lines is likely driven by radiation pressure on ions 
at UV resonance transitions. For the wind to be launched 
in the vicinity of the broad-line region, the quasar contin¬ 
uum must not only have signihcant power in the UV - the 
source of the line-driving photons - but not too much power 
in the X-ray. X-rays will strip atoms such as C and O of 
their electrons, thus dropping the absorption cross section 
dramatically. Overionized gas can therefore only be driven 
radiatively by Thomson scatt ering, which is much less e f- 
ficient than line drivi ng (e.g., lE'erett fc Ballantvnell2004l ~l. 
The elegant picture of [Murray et al.l 'i iggsl f that unifies the 
broad emission lines seen in all Type 1 quasars with the BAL 
wind directly observed in only ^ 20 % of optically selected 
quasars (e.g., [Hewett 8z Folt j [2003ll by generating broad 


© 2014 RAS 














2 Gallagher et al. 


emission lines in an equatorial disk wind is compelling, but 
has not been widely adopted in other models. In such mod¬ 
els, broad-line “clouds” are still invoked as the source of the 
broad-line emission, though the smoothness of the broad¬ 
line profiles and the requirement of some means of pressure- 
confinement to prevent the putative clouds from shredding 
rema in persistently unsolved problems i n the cloud picture 
(e.g.,[Murray fc Chiandl 19971 : though see lBottorff fc Ferlandl 
1200011 ■ 


An important part of active galactic nucleus (AGN) 
models that is even more uncertain is the so-called torus, 
a dusty molecular region presumably on the outskirts of the 
accretion disk. The inner radius of the torus is set by the dust 
sublimation radius, beyond which some fraction of the ac¬ 
cretion disk continuum is absorbed by grains and reradiated. 
Functionally, the torus must account for the ~ 30% of the 
integrated radiant quasar power that comes out in the near- 
to-mid-infrared. Furthermore, it must obscure the central 
continuum and broad-line region in a significant fraction of 
objects called Type 2 AGN. While a static toroidal structure 
(whether smooth or clump y, e.g., iKrolik fc Begelmanlll98^ : 
iNenkova et al.ll20n2l . 120081 1 beyond the dust sublimation ra¬ 
dius will serve these purposes, such a cold (100-1500 K) 
structure with the dust mass implied by the mid-infrared lu¬ 
minosity would quickly become gravitation ally unstable and 
collap se. This problem was addressed by iKonigl fc Karti j 
(Il994h with a magneto-hydrodynamically launched wind, 
shaped also by radiation p ressure on dust. Pr o posin g an 
alternate dynamical model, lElitzur fc ShlosmanI (l200fil l use 
magnetic fields to launch and confine dusty clouds from the 
plane of the accretion disk. ([Lawrence fc Elvig l2010l sug¬ 
gest an alternate picture where warped disks, originating 
from gas in the plane of the host galaxy, account for the 
large scale-height needed to obscure the central engine; the 
gravitational instabilities th at can give rise to su ch asymme¬ 
tries are further explored bv iHopk ins et al1l2012|. ) The dusty 
wind picture was expanded by [Keating et al.l ll2012l . here¬ 
after K12) by incorporating radiation driving on d ust into 
the m agneto-hydrodynamic (MHD) wind models of I Everett! 
1 I 2 OO 5 II . and then illuminating the wind with the central 
continuum to generate a model spectral energy distribution 
(SED). While the details of the shape of the output contin¬ 
uum still require some refining to match those observed, the 
overall power and general shape generated by the K12 ‘fidu¬ 
cial’ model is promising. This model uses as inputs the em¬ 
pirical SDSS composite quasar continuum (with the X-ray- 
to-optical flux ra tio, Qqx, set to -1.6 a s appropriate for the 
input luminosity: [Richards et ^|2006| L a black hole mass of 
Mbh = 10® Mq, an Eddington ratio of L/L^dd = 0.1, and a 
column density at the base of the wind of Nh,o = 10^® cm“^. 
We further assume a Milky W ay interstellar-medium dust 
distribution ([Praine fc Li[ [200^ 1. In this article, we expand 
on the results presented in K12 by exploring how key in¬ 
put model properties affect the structure of the wind and 
extracting other important observables that are generated 
by the wind model. Finally, we offer some predictions for 
observations that would further test this picture. 


2 THE DUSTY WIND MODEL 

Following [Blandford fc Pavnel ([l982[') and [Konigl fc Karti j 
|[l994ll . we model the torus as a self-similar dusty wind 
driven by MHD forces and radiation pressure. The MHD 
radiativ e wind code has been described in more detail pre¬ 
viously |[Everett[[2005[ : [Keating et al.[[^012[ l , and we give only 
a brief outline here. First, a purely MHD-driven wind solu¬ 
tion is calculated given a set of initial conditions input by 
the user. Next, clo udy (version 06.02.09b, as described by 
[Ferland et al.l[l998l l is called to calculate the radiation field 
and dust opacity in the wind, allowing for the model to de¬ 
termine the radiative acceleration of the wind due to the 
central accretion disk continuum. A new MHD solution is 
calculated for the wind, now taking into account the radia¬ 
tive acceleration, and the process is iterated until the wind 
converges to a stable geometry. 

Operationally, our assumption of self-similarity means 
that several important quantities, including magnetic field 
strength and mass density, scale with spherical radius; this 
simplifies the calculation of the solution to the MHD equa¬ 
tions considerably. Furthermore, it allows for the incorpora¬ 
tion of radiation pressure from the disk by simply reducing 
the effective gravitational potential. Importantly, it allows 
us to use one streamline as an approximate template for the 
entire wind, simply scaling appropriately inwards or out¬ 
wards to find properties of the specific region of the wind 
under consideration. In addition, these assumptions allow 
us to comparatively quickly estimate the large-scale prop¬ 
erties of these winds. The relevant input parameters and 
physical properties of the specific models discussed below 
are presented in Table [T] 


2.1 Density Distribution 


The final, converged MHD radiative wind model therefore 
provides the mass density of the wind as a function of r, the 
radial distance from the axis of rotation of the disk, and z, 
the vertical distance from the plane of the disk. In more de¬ 
tail, the model supplies the density along each streamline as 
a function of the mass density at the base of the streamline 
and of polar angle. The density of the wind at any arbitrary 
point is then relatively straightforward to calculate as a con¬ 
sequence of the self-similarity assumption. Because we know 
the shape of the wind streamlines, we can use the shape of a 
single streamline to trace back from an arbitrary point to the 
disk-launching radius for any wind streamline. From there, 
given a number density of particles on the disk at the inner 
launch radius and t he fact that density is assu med to be pro¬ 
portional to ® ([Blandford fc Pa^^[l982ll . we calculate 
the number density at the launch radius of that particular 
streamline. To convert from number density to mass density, 
we assign the proton mass to the average particle. 

In Figure [ 1 ] we present a 2D color map of the number 
density of one quadrant of the dusty wind as a function of r 
and z for the fiducial model. The inner wind launch radius 
is set by the dust sublimation radius. The notable thinness 
of the wind is a consequence of the model parameters. In 
general, the outer wind launch radius is given by: 


^ outer 



Nn,o \ ^ 


( 1 ) 
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Figure 1. A 2D cross-section of the number density map of the 
fiducial dusty wind model with A’h,o = 10^^ cm Mbh = 
10® 5 and L/LEdd = 0.1 (model 1 in Table The color 

indicates number density according to the scale bar on the right. 
The number density drops sharply as a function of height, z, as 
the wind accelerates. Dashed straight lines indicate the values of 
the optical depth r at HaA6563, as labeled. 


for user-specified values of the inner launch radius rin, the 
number density at the inner launch radius riin, and the col¬ 
umn density along the base of the wind A'’h,o- 


2.2 Defining the Covering Fraction 

Because the dust is entrained in the wind, E{B — V) (cal¬ 
culated explicitly by cloudy) is a smooth function of incli¬ 
nation angle, monotonically increasing toward the plane of 
the disk where the polar angle 9 — 90°. We see this contin¬ 
uous distribution of density as a strength of the dynamical 
model presented, however, these values offer no obvious in¬ 
clination angle to choose as the boundary for the wind, and 
the number density map does not show any obvious feature 
that would define a ‘height’ for the torus. The optical depth 
through the wind is also a smooth function of 6 for a wide 
range of wavelengths, and similarly provides no obvious an¬ 
gle to mark the boundary of the wind. We have therefore 
chosen the angle at which the dust reaches an optical depth 
of 5 for light at a wavelength of 6563 A (the wavelength of 
Ha emission) to be a conservative definition of the open¬ 
ing angle (i.e., the inclination angle up to which the central 
engine may be detectable with high signal-to-noise ratio 
spectra to an observer) of the dusty wind. At this value 
of r, only 0.7% of the broad Ha emission is transmitted 
through the wind, and so a quasar viewed from this inclina¬ 
tion angle (or closer to the disk) would be classified as an 
optical narrow-line object. For reference, this value of thc 
corresponds to E{B — F) rs 1.6 and Av ~ 5.5 (cloudy as¬ 
sumes Rv = 3.4). For our fiducial model, 0open = 69°; the 
opening angles for all of the models are listed in Table [T] 
A value of ^open = 69° gives a dusty wind covering fraction 
of cos(0open) =35%. As expected, this is approximately the 
fraction of the quasar bolometric luminosity that is emitted 
from 2-100/im in the predicted SEDs from K12. For com¬ 
pleteness, we also include the less restrictive opening angle 
values for thq, = 3 (corresponding to 5% transmission of the 
Ha line flux and Av ~ 3.3) which is more relevant for low 
signal-to-noise ratio optical spectra typical of extragalactic 
survey data. 



Figure 2. A cross section of the wind, as in Figure[T] comparing 
the profiles of several wind models. Shown here is the fiducial 
model (black, labelled “A” ; model 1 in Table [TJ, with A^h.q = 
10^® cm“^, T/TEdd = 0.1, and Mbh = 10® Mq; a similar model 
with Mbh = 5 x 10® Mq (red, “B”; model 2 in Table [TJ; and 
a third model with the same parameters as the fiducial model 
except that T/TEdd = 0.01 (blue, “C”; model 3 in Table [TJ. 
The dimensions of the model B streamlines have been scaled by 
l/v/5 to directly compare the shape of the outflow to the other 
models. Increasing the black hole mass by a factor of 5 (keeping 
L/LEdd constant) does not significantly change the shape of the 
wind, but reducing T/T/Edd by a factor of 10 moves the inner 
launch radius in by a factor of -v/TO (as expected), and causes the 
inner wind streamlines to be more vertical. This would have the 
effect of reducing the projected cross-section of the base of the 
inner wind (which is optically thick to near-IR emission; see Fig. 
7 in K12) compared to what is seen in the fiducial model. We 
propose this geometric explanation for the increased prominence 
of the 3-5 fim bump in luminous quasars. 


3 DEPENDENCE OF THE WIND ON INPUT 
PARAMETERS 

3.1 The Effect of T/TEdd 

Within the context of our model, we can explore variations 
in windy torus geometries across the AGN population by 
investigating the effects of changing luminosity, L /Z/Edd, and 
Mbh on the structure of the wind. Figure [2] illustrates the 
shape of the fiducial wind model (labeled as A; model 1 in 
Table [TJ compared to two other models. Model B (model 2 
in Table [TJ has a black hole mass of 5 x 10® Mq, five times 
that of the fiducial model. Because the value of L/L-Edd is 
fixed at 0.1, this model is 5 times more luminous than the 
fiducial model. The resulting shape of the wind is nearly 
identical to that of the fiducial model; the streamlines are 
slightly more vertical. Model C (model 3 in Table [TJ has a 
black hole mass of 10® Mq- the same as the fiducial model 
- but with L/L^dd dialed down to 0.01. The inner launch 
radius moves in by a factor of VTO, as expected with an 
scaling of rin, and the inner streamlines are more vertical as 
radiation pressure is relatively less important than at higher 
Eddington ratios. The slight change in wind structure from 
altering Mbh but keeping L/L^dd constant compared with 
the larger change from decreasing T/I/Edd by a factor of 10 
shows that the Eddington ratio is a primary driver given our 
modelled wind properties. 
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3.2 The Effect of fJ—field Strength 

The organized, poloidal magnetic fields in our fiducial MHD- 
wind model are reasonably strong, with mass-to-magnetic 
flux ratios of k = 0.03, implying considerable magnetic fields 
strengths of ~ 0.235 G at the innermost streamline. While 
an ordered magnetic field is required at the base of the wind 
to lift material into the line of sight of the accretion disk 
continuum, radiation pressure on its own may effectively ac¬ 
celerate material after that point. We explored the effect 
of dialing down the strength of the magnetic field. Because 
K oc -g, this means increasing k, and so we ran models with 
k= 0.04, 0.06, and 0.12 for values of LjL-Edd = 0.10 (models 
4-6 in Table [T|. Changing B to be up to 4x weaker had 
a significant effect on the distribution of E{B — V) values 
over the sky. We found that stronger fields (smaller values 
of k) result in higher E{B — V) values over a larger frac¬ 
tion of the sky; the Type 2 boundary for relatively stronger 
fields occurs at lower inclination angles ( 69° vs. 74° for 
K = 0.03 and 0.12, respectively). These results are a direct 
consequence of radiation driving dominating when the fields 
are weaker, thus leading to more radial outflows. 

We also ran the different field-strength models for lower 
T/Tfidd values of 0.05 (vs. 0.10 for the fiducial model), and 
found little effect on the cumulative E{B — V) distribution 
(see Fig. [3] see models 7-9 in Table [TJ. The lower L/L^dd 
values also lead to smaller rsub values (and thus rin values 
for the wind), but the E{B — V) distribution appears to 
respond most sensitively to increasing k because radiation 
pressure consequently becomes more dominant. At mass- 
to-magnetic flux ratios of k ~ 1 the MHD wind will fail, 
and without an alt erate source of vertical lift (e.g ., infrared 
radiation pressure; iDorodnitsvn fc Kallmanll2012l i radiative 
driving will also not be effective because material from the 
outer disk will not be exposed to the light from the inner 
disk. 


3.3 The Effect of X-ray Weakness 

Within the context of the dusty wind model, we are also 
interested in exploring the effect of the well-known anti¬ 
correlation between Oox (the 2 keV to 2500 A spectral index) 
and I/uv whereby more UV-luminous quasars emit a smaller 
fraction of their power in the X-rays. For our fiducial model, 
the value of Oox is —1.6; we have also run models assuming 
X-ray-loud (oox = —1.1) and X-ray weak (qox = —2.1) ver¬ 
sions of the illuminating SED. Figure 0] shows the opening 
angle, 0open, as a function of inner radius, rin. For this large 
range of Oox values (-2.1 to -1.1 encompasses the extremes 
of the observed distribution), the effect of X-ray loudness 
on the opening angle is not large, although two of the X- 
ray louder models do have slightly larger opening angles 
{9 ~ 75°) compared to the X-ray quiet ones {9 ^ 73). This 
result is reasonable, because unlike the ionization parameter 
of the associated gas, dust driving depends primarily on the 
integrated power in the optical through X-ray, and less on 
the spectral shape. 

The models with altered Oox presented in K12 were not 
particularly realistic, as the X-ray continuum was modeled 
simply as a flat (in units of power vs. frequency) continuum 
above 0.29 keV whose normalization was set by Qox. In an 
effort to create a more realistic continuum, we instead in- 
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Figure 3. The cumulative covering fraction of the wind {Q/2tt) 
vs. extinction, E{B — V)', larger r2/27r values at low extinction 
show winds that cover less of the sky with significant extinc¬ 
tion. The dotted vertical line marks Ay = 1, and the dashed 
vertical line marks the Type 1/Type 2 boundary as defined by 
THa = 5. The mass-to-magnetic flux ratios, k, for each model are 
as labelled; with thick (thin) curves representing L/L^dd = 0.1 
(0.05). The fiducial model curve (marked ‘fid’) has k = 0.03. The 
curves plotted are for models 1 and 4—9 in Table [T] Decreasing 
the strength of the magnetic field (increasing k) increases the 
importance of radiative driving and makes the wind more radial; 
consequently, the opening angle of the wind increases and the 
extinction to Type 1 objects decreases overall. 


elude a dust-free layer of shielding gas interior to the dusty 
wind. The shape of the ionizing continuum impinging on the 
dusty wind is therefore altered by gas absorption more re¬ 
alistically. To explore the effect of increasing the shielding 
gas thickness, we ran four additional versions of the fidu¬ 
cial model where only the column density of shielding gas 
has been altered with values of log(A^H, shield) = 19, 23, 24, 
and 25 cm“^ (models 10-13 in Table[TJ. The fiducial model 
has log(AH,shieid)=18 cm“^. When we calculate the cumu¬ 
lative distribution of E{B — V) as in Figure |3l we find that 
increasing the dust-free shield makes no difference in the 
shape of that curve, and the opening angles for the dusty 
wind remain the same as for the fiducial model. Therefore, 
the shielding gas does not strongly impact the covering frac¬ 
tion of the dusty-wind torus, consistent with the results from 
our X-ray weak continuum modeling. 


4 OBSERVATIONAL CONSEQUENCES OF 
THE DUSTY WIND MODEL 

The geometric structure, dust distribution, column density, 
and ionization state of material in the dusty wind model 
have direct observational consequences. In K12, we predicted 
mid-IR SEDs of the wind. Now, we are specifically interested 
in exploring if the dusty MHD wind can account for the 
ratio of Type 1 to Type 2 AGNs, the gross properties and 
prevalence of warm absorbers, and the near-IR bump from 
hot dust at ~ 1200 K in luminous AGNs. 
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Table 1. Properties of Dusty Wind Models. All models have Oox=-1.6 and Mbh = 10® Mq except for model 2 (fid_mbh8p7) with 
JWbh = 5 X 10^ M0. 


Model Name 


L/L^dd 

Mb 

H,shield 
log(cm“^) 

r9 

in 

(10^® cm) 

gd 

open 

(degrees) 
thq = 5 

^outflow 
{Mq yr-1) 

"^^term 

(10® km s“l) 

1. fiducial (A)® 

0.03 

0.10 

18 

4.57 

62.3 

69.4 

1.20 

2.94 

2. fid.mbh8p7 (B)® 

0.03 

0.10 

18 

10.6 

62.2 

69.3 

3.80 

5.00 

3. fid.edd001 (C)® 

0.03 

0.01 

18 

1.60 

60.4 

67.3 

0.99 

1.90 

4. sgl8_rsublg_k04 

0.04 

0.10 

18 

4.57 

63.7 

70.0 

1.21 

2.89 

5. sgl8_rsublg_k06 

0.06 

0.10 

18 

4.57 

66.9 

71.7 

1.21 

2.89 

6. sgl8_rsublg_kl2 

0.12 

0.10 

18 

4.57 

74.0 

76.5 

1.22 

2.77 

7. sgl8_rsubsm_k04_edd05 

0.04 

0.05 

18 

3.23 

63.0 

69.1 

1.07 

2.81 

8. sgl8_rsubsm_k06_edd05 

0.06 

0.05 

18 

3.23 

66.5 

71.2 

1.07 

2.78 

9. sgl8_rsubsm_kl2_edd05 

0.12 

0.05 

18 

3.23 

73.8 

76.4 

1.08 

2.77 

10. sgl9.rsub4pt570 

0.03 

0.10 

19 

4.57 

62.3 

69.4 

1.20 

2.94 

11. sg23j-sub4pt480 

0.03 

0.10 

23 

4.48 

62.3 

69.4 

1.19 

2.95 

12. sg24_rsub4pt320 

0.03 

0.10 

24 

4.32 

62.3 

69.4 

1.18 

2.97 

13. sg25_rsub3pt900 

0.03 

0.10 

25 

3.90 

62.3 

69.4 

1.13 

3.03 


“Input mass-to-magnetic flux ratio. A smaller value indicates a stronger poloidal S-field. *'Input column density of the 
shielding (dust-free) gas at the base of the wind. “The dust sublimation radius at the base of the wind calculated based 
on the input parameters. This defines the launching radius of the innermost wind streamline. “^The opening angle of the 
dusty wind as defined for r = 3 and 5 at HaA6563. “The mass-outflow rate of the dusty wind. ^The terminal velocity of 
the dusty wind. ® Letters in parentheses refer to the model labels in Figure [2] 


4.1 Broad vs. Narrow Line AGN in the Dnsty 
Wind Picture 

Given the requirement that a successful torus model ob¬ 
scure the accretion-disk continuum and broad-line region 
emission, comparing our predicted fraction of Type 1 (op¬ 
tical broad-line) and Type 2 (optical narrow-line) quasars 
to observations is an important test of the models. In our 
models, the range of opening angles is relatively tight, and 
shows little dependence on model input parameters with the 
exception of the column density at the base of the wind (only 
models with Ah.o = 10^® cm“^ are shown in Fig. 0, and the 
strength of the magnetic field (see Figs.[3]and0. Specifically 
for the suite of models run by K12, typical opening angles 
using the thq = 5 criterion are Oopen = 69-75° (see Fig. 0; 
this corresponds to dusty-wind covering fractions (Type 2 
fractions) of = cos (Sopen)=26-36% or Type 1 to Type 2 
ratios of 1.8 to 2.8. In their su rvey of low - redshi ft, emission¬ 
line galaxies from the SDSS, iHao et ^ l|2005h found that 
Type 1 objects outnumber Type 2 objects by a factor of 2- 
4 at high luminosities. A similar study by ISimpsoril il2005l ) 
found Type 1 fractions of ~ 60% (Type 1/Type 2 = 1.5) at 
the luminous end (in terms of [O ill] luminosity). Our num¬ 
bers are consistent with (though narrower than) the range 
of 1.5 to 4 in these population studies. 

As demographic studies can be plagued by selection ef¬ 
fects, an alternate means of constraining the covering frac¬ 
tion of the dusty “torus” is to systematically measure the 
power in the IR-region of Type 1 quasar SEDs. This as¬ 
sumes that the dusty medium acts effectively as an efficient 
bolometer whose emitted power is simply the fraction of 
the sky covered by the dusty medium from the point of 
view of the accretion disk. While radiative transfer mod¬ 
eling shows this is not strictly true (cf. K12), it is a reason- 
ab le approximation . Alon g these lines, the recent analysis 
by iRoseboom et al.l (l2013l ) using WISE observations of lu¬ 
minous (~ 10^® erg s“^) Type 1 SDSS quasar determined 
that the covering fractions could be described by a log nor¬ 


mal distribution with mean of 391^4%. This correponds to 
an intrinsic Type 1/Type 2 ratio ranging from 0.61 to 3.0; 
the range from our suite of models of 1.8 to 2.8 is again 
within (though narrower than) the span of empirical values. 

A Type 2 fraction that decreases with increasing lu¬ 
minosity does not follow simply from our model in the 
parameter range we have explored, as the wind covering 
fraction is not sensitive to only luminosity, but a combi¬ 
nation of L/I/Bdd and the mass-to-flux ratio of the mag¬ 
netic field. However, it should be kept in mind that there 
can be other explanations for a luminosity dependence for 
obscuration. For example, the host galaxy can also serve 
as a signi ficant source of ob s curation of the bro ad-line re¬ 
gion (e.g.. iRigbv et al.]|2006l : [juneau et al.ll201ll) . If accre¬ 
tion disks are randomly oriented with respect to their host 
galaxies as observed, the observed Type 2 fractions are in¬ 
flated with respect to the torus-covering fraction because of 
line-of-sight obscuration occurring from optic ally thick dust 
within the host galaxy (iKinnev et al.l I 2 OOOI) . Because the 
size of the BLR scales w ith luminosity (oc Z/° ®: lKasDi et all 
I 2 OOOI: iBentz et al.ll2009l) . more luminous broad line regions 
are larger and therefore harder to hide, even before invok¬ 
ing any sort of feedback mechanism that might clear the 
host galaxy ISM from the vicinity of the quasar. Further¬ 
more, weak broad lines from low-luminosity objects can 
be “hidden” by a stron g stellar continuum dMoran et al.l 
I 2 OO 2 I : [Martini et al.ll200d) : detecting these requires high S/N 
spectra and careful modeling and subtracting of the stellar 
galaxy emission (e.g.. lHo et al.lll9^ : iHao et al.ll2005l) . 

Here we are defining Type 1/Type 2 as broad- 
line/narrow-line with respect to Ha; however, the X-ray 
community often uses the detection of intrinsic X-ray ab¬ 
sorption with Ah > 10^^ cm“^ to distinguish obscured vs. 
unobscured quasars. The two criteria overlap, but do not de¬ 
fine identical populations, as dust extinction and gas absorp¬ 
tion do not always go hand-in-hand. The Ah value at what 
we have defined as our Type 1/2 boundary is ~ 10^^ cm“^. 
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(1 o’® cm) 


Figure 4. Top: The opening angle, ^open (degrees), as a func¬ 
tion of the wind launch radius, rin (10^® cm). All models have 
base column densities of A^h,o = 10^^ cm~^, the value of A^h,o 
that generates appro priate near-to-mid-infrared power in the SED 
llKeating et al.|[2012l) . The crosses, diamonds and asterisks indi¬ 
cate values of Qox as labeled in the legend. The diamonds are 
scaled to represent their L/LEdd values; the size shown in the 
legend is L/LEdd = 0.05. The models are those presented in K12; 
the bulk of them give opening angles of 69°—75°. The outlier with 
Qfox = —1.6 at 0open ~ 86° is the model with rin set artificially 
low to mimic the effect of increasing the dust sublimation tem¬ 
perature. Bottom: Same as panel above for the new models with 
variable mass-to-fiux ratios and variable shielding gas column 
densities. Diamonds are scaled to represent k. values from 0.03 
(smallest, i.e., strongest magnetic field strength) to 0.12 (largest; 
weakest magnetic field strength). Asterisks show k, = 0.03 models 
with a range of shield values (log(A^H shield) [cm~^]) as labeled; 
all other models have log(A^H, shield) = 18 cm~^. While decreas¬ 
ing the strength of the magnetic field increases 0open, increasing 
the column density of the shielding gas has little effect on ^open, 
although the launching radius decreases with increasing A^h shield- 


and therefore consistent with both the X-ray and optical 
definitions. In the wind model, this also sets a floor for the 
amount of intrinsic absorption that would be found in opti¬ 
cal Type 2 objects where the dusty wind is obscuring the 
broad line region, consistent w it h the Seyfert 2 X-ray ab¬ 
sorption study of iRisaliti et al. (Il999ll and observa tions of 
X-ray occultation events by Markowitz et al. I (EH). 


4.2 Ionic Absorption Through the Wind 


4-2.1 Warm Absorbers 

As seen in Figure[Tl wind density drops quickly as a function 
of height above the plane of the accretion disk. Because the 
wind is relatively collimated and narrow radially, this means 
that the line-of-sight column density also drops rapidly with 
decreasing inclination angle. The wind at moderate inclina¬ 
tion angles is therefore exposed to the high-energy contin¬ 
uum from the accretion disk and becomes highly ionized. 
Such highly ionized gas along the line-of-sight to the X-ray 
continuum is observed in m any AGN (~ 50% of Seyfert 
galaxies; e. g.. lRevnoldslfl997ll : the strongest atomic features 
of these so-called “warm absorbers” are the bound-free edges 
for O VII and O vill and the helium-like and hydrogen-like 
inner shell bound-bound transitions of these ions. 

Studies of the variability of X-ray absorption features 
have led to loose constraints on the location of warm ab¬ 
sorbers that are la rgely consistent with the torus location 
(e.g., within ~ 6 Pc: lKrongold et ^l2005l l. The specific iden¬ 
tification of the warm absorber with a hot phase within the 
vicinity of a molecular torus has been explor ed previously. 
Indeed, this was suggested bv iRevnoldsI lll997ll when ASCA 
observations revealed the prevalence of warm absorbers in 
Seyfert galaxies with some preferen ce for those with heav - 
ier optical extinction. Subsequently, iKrolik fc Kris j ll200ll ~) 
explored a multi-temperature thermal-wind origin for the 
warm absorber gas: they postulated that warm absorbers 
were evaporated off of the inner surface of the dusty torus. 

More recently, iDorodnitsvn fc KallmanI (l2012l 'l pre¬ 
sented a radiation-hydrodynamic model that incorporated 
infrared radiation pressure on dust to drive a slow wind (10- 
200 km s“^) off of the accretion disk beyond the dust subli¬ 
mation radius. The inner (from the point of view of the ac¬ 
cretion disk) layer of this slow, cold wind is illuminated by a 
bright X-ray continuum to form a hot, fast (400-800 km s“^) 
wind. In the parameter space they considered, UV radiation 
pressure on dust grains was not included, and the X-ray lu¬ 
minosity needed to power the wind was extremely high, up 
to I/x = 0.3I/Bdd, much larger than seen in luminous quasars 
(where Lx ~ O.OlLsdd is more typical). The opening angle 
for the Thomson-thick part of the torus was 72° — 75°, which 
is considerably smaller than our values (close to 90°). 

With high S/N X-ray gratings spectra, O vil and O vill 
lines and edges in absorption can be used to constrain the 
column densities of the X-ray absorbing oxygen ions. For 
reference, two high-quality examples of this type of analysis 
are based on the 9 00 ks Chandra observation of NGC 3783 
dKaspi et al.l 120021') and the simultaneous HST / Chandr a 
study of Mrk 279 (IGostantini et ^ ] |2007l: lArav et al.l liH) ; 
both Seyferts are X-ray bright with particularly strong warm 
absorber features. The best-fitting column densities from 
these studies range from 1-5 x 10^^ cm~^ (O vil) and 0.3- 
39 X lO’^^ (O vill). The more recent X-ray study of t he red¬ 
dened quasar IRAS 13349+2438 bv iLee et al.l ( 2013ll found 
similarly high column densities of 10^^“^® cm~^ for both 
ions. From the wealth of features in the HETGS spectrum, 
the authors favoured a model with an absorber with a con¬ 
tinuous distribution of ionization parameters, and suggested 
that it could b e an ionized atmosph ere above the putative 
torus as in the lKrolik fc Kri^ (l200ll 'l picture. 

The column densities from these detailed case studies 


© 2014 RAS, MNRAS 000.IT1ITT1 














































The Windy Torus 7 


of the strongest (and therefore - by selection - not typi¬ 
cal) warm absorbers are about an order of magnitude higher 
than we find in our fiducial model for O vill, and about a 
factor of 100 larger that our values for O vii (see Fig. [5]). 
Our X-ray weak modification to the fiducial model (with 
Oox = —2.1) however, gives values of O vill and vil that 
are both ~ 10^^ cm“^ over a large range of inclination an¬ 
gles. The X-ray loud (oox = —1-1) fiducial model predicts 
much lower values for the O vill and O vil column densities 
because most of the oxygen nuclei are completely stripped 
of electrons. While the distributions for the fiducial models 
with the X-ray modifications are all quite flat as a func¬ 
tion of inclination angle for 0 < 6^open, the k = 0.12 model 
(with the weakest B-field), has a much smaller range of an¬ 
gles over which the “warm absorber” gas would be detected 
(from 0 = 55° — 76°), but the hydrogen-like and helium¬ 
like oxygen column densities rise steeply as the inclination 
angle approaches 0open, and near that boundary have values 
consistent with the empirical values measured from the ex¬ 
treme X-ray spectra cited above (see Fig. El). Adding a layer 
of dust-free shielding gas up to Ah, shield = 10^® cm“^ to 
the fiducial model at most approximately doubles the O vii 
and O VIII column densities. Given the range of values ob¬ 
tained from the models, we consider the high ionization gas 
within the MHD-|-radiation-driven wind a plausible expla¬ 
nation for more typical (lower column density) X-ray warm 
absorbers. All of our models by design have “dusty warm 
absorbers”; the dust is collisionally coupled to the gas which 
is launched and collimated by magnetic fields to large scale 
heights where it can be exposed to the X-ray continuum and 
cover a large fraction of the sky. The model in its present 
form checks (via cloudy) that the gas temperature is con¬ 
sistent with dust survival along the entirety of a streamline 
to determine the innermost boundary of the wind. Nonethe¬ 
less, incorporating the data from Figs. El and EJi, there is a 
large range of inclination angles {0 = 15-45°) in the fidu¬ 
cial model with a substantial warm absorber column density 
of Aoviii ~ 10^® cm“^ without severe reddening {Ay < 1). 
The models with weaker B-fields (see Fig.Et*) generate dusty 
warm absorbers exclusively. 

It is not obvious how one could distinguish observation- 
ally between the different proposed scenarios for the warm 
absorber origin with present data, but future missions such 
as Athena+ that will better detect and resolve the velocity 
structure of X-ray absorption f eatures (including dust sig¬ 
natures; e.g., iLee fc Ravei2005h may shed some light on the 
relationship of the warm absorbers to a dusty wind. Our 
models currently provide only integrated column densities 
of lines-of-sight through the dusty wind, and so we cannot 
directly compare absorption profiles as a function of outflow 
velocity with observed spectra in the UV and X-ray. To pro¬ 
vide some reference, the maximum vertical outflow velocities 
are a few 10® km s“® (see Table [TJ. Warm absorbers typ¬ 
ically have multiple velocity components present that can 
vary independently, and the radiatively driven dusty wind 
could plausibly account for some of them. However, explain¬ 
ing the details of individual X-ray spectra frequently re¬ 
quires some sort of inhomogeneous structure, and we do not 
claim that the MHD-wind in this iteration of our model is 
sufficient to account for all of the spectral features. Nonethe¬ 
less, the presence of highly ionized gas in the dusty outflow 


occurs naturally within the MHD plus UV/optical and X- 
ray radiation-pressure driven model. 


4-2.2 Low Ionization Gas 


If we use Mg ii as a tracer of low ionization gas, it is apparent 
from Figure E] that significant column densities of this ion 
are not visible in Type 1 objects; the column density of Mg ll 
only rises steeply at inclination angles larger than the open¬ 
ing angle. In this case, the dusty wind alone cannot give rise 
to the low ionization gas responsible for broad Mg ll absorp¬ 
tion lines seen in LoBAL quasars. However, in the fiducial 
model, we have not included a layer of dust-free shielding gas 
or the UV BAL outflow, as would of course be required for 
a robust model of broad absorption line quasars. Along the 
line-of-sight from the X-ray continuum from the inner ac¬ 
cretion disk, BAL quasars have column densities of at least 
10^^ cm“^ and low-i onization BAL quasar s may be either 
Compton-thick ( e.g., Gnllagher et al.fl2006h or intrinsically 
X-ray weak ('e.g.. lLuo et al.l 2013ll with observed Qox values 
of — 2.0. (Our addition of MHD-driven shielding gas for 
some model runs does not include UV radiation-driving, and 
so it is not a realistic model for a BAL quasar.) 

Relatively rare among optically selected quasars, 
LoBAL quasars h ave dust-reddened continua that are intrin¬ 
sicall y quite blue (ISpravberrv fc Foltj|l993 : iReichard et al.l 
2003li. and may b e evolutionar i ly yo ung systems (e.g.. 


Farrah et al.l 120051 : lUrrutia et al.l l2009l L Alternatively, all 


BAL quasars may have a LoBAL region that is only probed 
through some fraction of the lines of sight through the wind. 
This is expected in a radiatively driven wind launched from 
an accretion disk given the stratified structure of the broad 
line region where higher ionization gas (probed by C iv) 
originates at smaller radii than lower ionization gas (probed 
by Mg ll). The similarity of the high ionization and low- 
ionization BAL velocity structure (at least at low velocities; 
IVoit et ^Il993li is consistent with both features being found 
within the same structure launched from a range of radii 
close to the source of the UV-optical continuum emission 
and well within the dust sublimation radius. The dust ex¬ 
tinction which is common in LoBAL quasars may however 
arise from a dust-driven wind as presented in this paper. 
The concurrence of extinction with LoBALs in a single ob¬ 
ject therefore would not indicate that there is necessarily 
dust entrained in the LoBAL outflow; the large inclination 
angles required to pass through the LoBAL gas are also con¬ 
sistent with larger values of E{B — V) from the dusty wind. 


4.3 Hot Grains in Dusty Winds 


Within the context of the wind model, we would also 
like to investigate the empirical observation that the 
strength of the 3-5 /rm bump - the manifestation of 
a larger contribution to the near-to-mid-infrared SED 
from hot dust - is seen to increase with increas¬ 
ing quasar lumino s ity (e.g., |Edelson__&Wlalkan| [lOsB; 
Richards_et_^ 200^ ^aUagher et al.ll2007l : Krawczvk et al.l 


20131 : Roseboom et al.l 2013l L In~ our picture, the inner 

streamline carries the hottest dust, and it becomes more 
visible as L/I/Edd increases because the relative importance 
of radiative driving goes up and the wind is pushed closer to 
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Inclination Angle [degrees] 


Figure 5. Top (a): Column density vs. inclination angle of abun¬ 
dant ions characteristic of both high ionization (e.g., O vii and 
O viii) and low-ionization gas (e.g., H i and Mg ii) for the fiducial 
K12 model. Bottom (b): The same plot as above for the lowest 
B-field model with k. = 0.12. Because of the increased importance 
of radiative driving, the ionized gas covers a significantly smaller 
fraction of the sky in this model, and the O vii and O viii column 
densities increase more steeply and to larger values than in the 
fiducial model with k. = 0.03. In both panels, the column densities 
of high-ionization gas are largely consistent with constraints from 
warm absorber studies. The low-ionization species have negligible 
column densities at inclination angles smaller than the opening 
angle (marked with a vertical dashed line). 

the equatorial plane ([Gallagher et al.ll2007l . K12). We pre¬ 
dict that this would have the effect of increasing the promi¬ 
nence of the 3-5 /rm emission, as the projected cross section 
of the surface of the inner streamline from a typical Type 1 
orientation angle will increase. The dusty wind becomes op¬ 
tically thick to near-infrared emission at large inclination 
angles (see Fig. 7 in K12). Consequently, within the dusty 
wind paradigm, we expect that the claimed correlation of the 
strength of the 3-5 /rm emission with luminosity is in fact 
a correlation with Eddington ratio. (For luminous quasars, 
luminosity and L /I/Edd are of course correlated given the re¬ 
stricted range of obs erved black hole masses in f lux-limited 
quasar surveys; e.g., IVestergaard k. Osmerllioogl .'l 

A second means of boosting the 3-5 fj,m emission is to 
decrease the inner launch radius in the hducial model by 
about one-third by increasing the dust sublimation temper¬ 
ature to approximate the effect of a graphite-driven wind 
(K12). This decrease in the inner launch radius causes the 
wind to be notably broader in our model similar to what we 
saw with the decrease in L/Lfidd and as expected from Equa¬ 
tion [T] However, the shape of the streamlines is very close 


to that of Models A and B. Because the bulk of the near- 
infrared emission is coming from the innermost streamlines, 
it therefore follows that a smaller rin for a given L /LEdd gen¬ 
erates more near-IR emission because of the higher dust tem¬ 
peratures (K12). However, the opening angle for the wind 
becomes quite large, 0open = 86° (see Fig. U)). If this is the 
cause of the luminosity dependence of the 3-5 jam emis¬ 
sion, this implies some mechanism for decreasing the dusty 
wind launching radius as a function of luminosity which is 
challenging to arrange. We therefore prefer the geometric 
explanation above. 

4.4 Images of the Wind 

As luminous quasars are typically farther away and there¬ 
fore fainter and of smaller angular size than nearby Seyferts, 
they are beyond the reach of the current generation of in¬ 
terferometers. Therefore we can only compare our simulated 
model structure to observations of the nearest Seyfert galax¬ 
ies, which are by nature considerably less luminous (because 
of lower Mbh masses or lower Eddington ratios or both). 
A simulated high-resolution image of the wind at 9.5 jim 
heated by the central continuum is shown in Fig. The 
hourglass shape is characteristic of MHD winds, and rem¬ 
iniscent of the outflows seen in young stars. It illustrates 
clearly that the base of the wind and the inner stream¬ 
lines dominate the emission at this wavelength; the bright¬ 
est 9.5 jam emission regions would correspond to the largest 
extinction values in the optical. Int riguingly, 8 t o 13 jam 
interferometric observations have led iHonig et al.l ||2012|') to 
infer a similar structure and characteristic size scale (~ 2 pc) 
for the dust distribution in the Seyfert 2 galaxy NGC 424. 
The largest compilation of mid-infrared interferometry to 
date shows that a large fraction of this emission is unre¬ 
solved, or at scales of < 0.1 to < 10 pc fo r a sample of 23 
local Seyfert nuclei dBurtscher et al.ll2013li . Reverberation¬ 
mapping using the response of the near-IR continuum to 
changes in the optical continuum of a several Seyfert galax¬ 
ies has placed the inner ‘torus’ radius beyond the broad-line 
region and consistent with th e expected dust-sublimation 
radius l|Kishimoto et all 120131 '). However, more direct com¬ 
parisons of the dusty wind models with observations await 
the light-collecting power and angular resolution (~ 0.001” 
at 2^m) expected in the era of 30-m telescopes. 

4.5 Constraints on Organized Magnetic Fields 

The success of the dusty wind model that we have pre¬ 
sented is predicated on the existence of relatively strong, 
organized magnetic fields that thread the outer accretion 
disk of quasars. Certainly, radio-loud AGNs have such fields 
in the inner disk; radio synchrotron emission from jets 
requires them. The picture for radio-quiet AGNs is less 
clear. There is some observational evidence for ordered 
magnetic fields within a few parsecs of the central en¬ 
gine ofJocal_AGNs_as_de^cted_jn near-infrared polarime- 
try. iLopez-Rodriguez et al.l ll2013h measured polarization 
of a few percent in the band for the Type 2 AGN 
IC 5063. They interpreted the polarization as being due 
to dichroic absorption of dust grains aligned by magnetic 
fields in the nuclear region. The polarization properties 
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Figure 6. A simulated high-r esolution 9 .5 fim image of the fidu¬ 
cial model generated by mc3d The hot “disk”-shaped 

structure is approximately 2 pc in radius. The bulk of the emis¬ 
sion comes from the base of the wind and the inner streamlines 
heated by the accretion disk continuum thus giving the MHD 
wind its characteristic hour-glass shape. 


thus constrain the B—field in the plane of the sky to be 
10-130 mG. However, their angular resolution is signifi¬ 
cantly larger (~ 260 pc) than the torus region, and thus 
the field may be associated with the interstellar medium 
of the host galaxy near the nucleus, rather than the torus 
itself. Similarly, from an analysis of Zeeman splitting in 
OH megamaser emissi on in ultra-luminous infrared galaxies, 
iR.obishaw et al.l l|2008l l estimate magnetic fields strengths of 
mG are present throughout the interstellar medium in these 
starbursting galaxies. On smaller scales (a few tenths of a 
parsec), the lack of detections of circular or linear polar¬ 
ization in measurements of the H 2 O molecular masers in 
Type 2 AGNs gives upper limits of 10-100 mG for the 
magnetic fields in the masing region (| Modiaz et al.l l200fil : 
IVlemmings et aP l2007l : iMcCalluni et a.1.1 2007fl . while our 
fiducial model has a field strength at rin near the disk sur¬ 
face of ~ 0.2G. B-field strengths in our k = 0.12 model are 
a factor of four smaller at 50 mG near the disk surface, and 
therefore still viable given these empirical constraints. Fur¬ 
thermore, the low luminosity Type 2 AGNs being studied 
(NGC 4258, NGG 3079, and the Gircinus galaxy) may have 
central regions that are quite different from Type 1 quasars, 
and so additional empirical constraints on poloidal B-field 
strengths within the inner parsecs of quasars would be valu¬ 
able. 


5 SUMMARY AND CONCLUSIONS 

We explore the properties of the wind models first presented 
in K12, and extract key parameters including the covering 


fraction of the optically thick wind, the column densities as 
a function of opening angle of several ions, and the distribu¬ 
tion of extinction values for Type 1 quasars. In addition, we 
explore the effect of decreasing the strength of the poloidal 
magnetic field and increasing the column density of the dust- 
free shielding gas interior to the dusty wind. Altering key in¬ 
put parameters of the model allows us to explore the effects 
of Mbh, L/L^dd, the shape of the ionizing continuum, the 
mass-to-magnetic flux ratio, and the column density of the 
shielding gas on the wind. We find the following empirical 
consequences of this model: 


• Using our definition of r < 5 at Ha to define broad-line 
quasars, the dusty wind model predicts a Type 1 fraction of 
~ 65% for our fiducial model (Table [T] and Fig. HI. Type 2 
quasars in our model have line-of-sight column densities of 
Ah ~ 10 ^^ cm“^ (Fig-Ejl. Both of these values are roughly 
consistent with observed Type 1 fractions and line-of-sight 
column densities found in optical and X-ray surveys of lu¬ 
minous AGNs and from mid-IR SED analyses. 

• While the presence of a UV line-driven wind as seen 
in broad absorption line quasars apparently requires a blue 
UV continuum and weak X-ray emission (either as a result 
of absorption or intrinsic weakness), the dust-driven wind is 
much less sensitive to the shape of the continuum. Mid-IR 
continuum emission might therefore be quite similar among 
classes of quasars with a range of ionizing SED shapes, and 
consequently a more accurate measure of bolometric lumi¬ 
nosity. 

• We present predictions for the cumulative distribution 
of E{B — V) values for quasars (Fig.|3}- The normalization of 
these curves depends primarily on k, the mass-to-magnetic 
flux ratio, with higher values of k (weaker B-fields) resulting 
in slightly larger opening angles, and significantly less dust 
extinction over a wide range of opening angles. This is a 
consequence of radiation pressure dominating over magnetic 
driving which results in more radial outflows. The effect of 
changing k has a larger effect on the E{B — V) distribution 
than changing L/I/Edd within the range we considered. 

• The column density through the wind drops rapidly 

with height from the disk causing the ionization parame¬ 
ter of the wind to increase dramatically. As a consequence, 
high ionization atomic species are present at a large range 
of inclination angles. The modeled values for the column 
densities of O vill through the wind are ~ ]^q 16-18 jjg_ 

pending on the properties of the ionizing continuum and the 
value of K (Eig.[SJ. This range is low for the strongest X-ray 
warm absorbers in the best-studied local Seyfert galaxies, 
but consistent with more typical warm absorbers. Improved 
measurements of the profiles of warm absorber features with 
next generation X-ray observatories may determine if the 
dusty wind model presented here remains a viable explana¬ 
tion for X-ray warm absorbers. 

• The hot dust emission present in the near-IR SEDs of 
luminous quasars could be associated with the innermost 
surface of the dusty wind, where grains with higher subli¬ 
mation temperatures than silicates can survive. We predict 
that the strength of this emission feature will correlate with 
L/LEdd as higher L/I/Edd flows are more radial, and the 
cross-section of the inner streamline visible to the observer 
at typical inclination angles becomes larger. Scatter could be 
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introduced in this correlation by a distribution of magnetic 
field strengths at given values of L/L^dd- 

• MHD+radiation wind driven models of the AGN ‘torus’ 
predict hour-glass shapes for the mid-IR emission which may 
be routinely detectable in nearby AGN with improvements 
in near-to-mid IR interferometry and the next generation of 
30-m-class telescopes. The shape of the innermost streamline 
(which generates the bulk of the emission) is most sensitive 
to the mass-to-magnetic flux ratio, and secondly to L/L-Edd- 

Getting the dusty-wind model to the point where it can 
plausibly account for the detailed structure in near-to-mid- 
infrared quasar SEDs requires additional refinement. In par¬ 
ticular, the current single-zone dust structure is not realistic, 
and the dust composition and grain-size distribution likely 
varies from that of the Milky Way ISM as assumed in our 
models. The need for different dust properties for quasar 
environs is already implied by observations indicating that 
the shape of the extinction curve t o most quasars is mor e 
consistent with SMG-type dust ('e.g.. lKrawczvk et al.ll2014h . 

While clouds are not likely to persist in the dy- 
namic environment arou nd a supermassive black hole (e.g., 
ISchartmann et ^l201lll . it is also unlikely that any wind 
structure is completely smooth. Including inhomogeneities 
in the wind is therefore an obvious next step, in particular 
to more realistically predict ionic column densities. 

More generally, the viability of the dusty MHD disk 
wind paradigm requires an investigation into whether rela¬ 
tively strong ( Xi 50 mG), large-scale poloidal magnetic fields 
persist within the central parsecs of quasar host galaxies. 
Certainly quasars with jets have large-scale magnetic fields; 
radio synchrotron emission is observed. In the near future, 
the available sensitivity and spatial resolution of ALMA and 
the Jansky VLA are promising for placing meaningful con¬ 
straints on the strengths of organized fields within the cen¬ 
tral few parsecs of radio-quiet active galactic nuclei. 
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